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Abstract 

The decay of 2-gluon colour singlets in quarks: 2g ^ qq + 2q2q has been sim- 
ulated with the Monte-Carlo method, taking into account an effective 1-gluon ex- 
change interaction between the emitted quarks, which was folded with a 2-gluon 
density determined self-consistently. 2-gluon densities were found with different 
radii, which correspond to 0+^ glueballs of the size of light qq, ss, cc, bb and heavier 
qq systems. Binding potentials between the two gluons have been deduced, which 
are consistent with the confinement potential from lattice results. However, self- 
consistency for the deduction of 2-gluon densities requires massless (or very light) 
quarks for all flavours. The masses are given by the binding energies of quarks 
and gluons, yielding excitation spectra of glueballs and J/^* and T states 
consistent with observation. The sum of q-q potentials yields a strong coupling Og 
consistent with the available data up to large momenta. 

The nucleon is described by a gluonium state coupled to 3 valence quarks, yielding 
ground state and radial excitations consistent with experiment. Finally, we discuss 
the compressibility of the nucleon and relate it to that of nuclear matter. 

PACS numbers: 12.38.Aw, 12.39.Mk, 14.20.Dh, 14.40.-n 

1. Introduction 

The two key problems in the understanding of the strong interaction are the confinement of 
quarks and gluons and the origin of mass, both related to the non-perturbative structure 
of quantum chromodynamics (QCD). A linearly rising confinement potential between 
quarks has been derived in potential models [H [2] and lattice QCD simulations [3], H], 
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but its origin is not well understood. The mass term in the QCD Lagrangian is also not 
understood, but for the generation of mass a coupling of the quarks to a scalar Higgs 
background field has been proposed. Finite quark masses give rise to the axion problem, 
which has not been solved. 

For the description of QCD in the non-perturbative regime mainly two non-perturbative 
methods have been applied, solutions of Dyson-Schwinger equations [5j and lattice QCD |6j, 
which solves the QCD equations by path integral methods on a space-time lattice. 

2. Deduction of 2-gluon densities 

In this paper a new phenomenological method is presented, which starts from the con- 
jecture, that the non-Abelian structure of QCD may generate bound 2-gluon systems, 
which decay into qq pairs. For the description of such bound states $ we write the radial 
wave functions in the form ?/'$(r = fi — = [ipi{ri) 1^2(^2)], where V'i(^) ^-^'g the radial 
wave functions of the two gluons. To investigate the properties of such 2-gluon systems 
we studied the decay 2g ^ qq + 2q2q with an attractive interaction between the emitted 
quarks. 

Assuming an effective 1-gluon exchange interaction Vig{R) = —ag/R between the emitted 
quarks with relative distance i? = jr'j — fj\, the decay from a 2-gluon system 2g —* (qq)"" 
requires a modification of the free q-q interaction by the density of the 2-gluon system, 
which may be expressed by a folding integral 

V,,{R) = I dfp^if) V,g{R-r) , (1) 

where p<^{f) is the 2-gluon density p<^{f) = |'?/'<j,(r )p. 

It is interesting to note, that for a spherical density the Fourier transform of eq. ([1]) to 
momentum (Q) space yields 

VM) = pUQ) , (2) 

where pq>{Q) = Att J r^dr jo{Qr) p<i>{r). Comparing this with the standard 1-gluon ex- 
change force yields a Q-dependent strong coupling as{Q) = as p<i>{Q), which is qualita- 
tively consistent with the known fact of a "running" of as{Q) and the condition as(Q)— >■ 
for Q^oo (asymptotic freedom). 
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Further, finite size of tlie decaying 2-gluon system liave been taken into account, as well 
as the fact, that the decay into 2g2g favours relative angular momentum L=0 between the 
emitted quarks, whereas for the decay in qq the outgoing quarks are in a relative p-state 
(L=l). Details are given in ref. [7]. By relativistic Fourier transformation [8] the effective 
interaction ([1]) can be transformed to momentum space 

V,,{Q') =471 J RHR UQ'R) V,g{R) , (3) 
with Q' = Q Jl + [Q^/4m|] and m$ being the mass of the 2-gluon system. 




Figure 1: Upper part: Resulting 2-gluon momentum distributions (multiplied by Q^) for 
decay in qq and 2g2g and sum. Lower part: Deduced 2-gluon density with estimated error 
band. 

Monte-Carlo simulations of gluon-gluon scattering have been performed in fully relativistic 
kinematics, in which the 2 gluons in the final state can decay in qq and 2g2g (using 
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massless quarks). The potential Vqq{Ap) ([3]) has been used as a weight function between 
the outgoing quarks (with the relative momenta Ap = Pi—pj). Resulting gluon momentum 
distributions dqq{Q) and d2q2q{Q) for decay into qq and 2g2g were generated. Their sum 
Dq>{Q') = dqq{Q') + d2q2q{Q') cau bc related to the radial density p$(r) of the 2-gluon 
system 



with Q' as in eq. ([3]). 

The condition, that p$(r) in the interaction ([3]) and in eq. (jlj) should be the same, allowed 
us to determine this density. Resulting momentum distributions dqq{Q) and d2q2q{Q) for 
a self-consistent solution with < >^0.5 fm^ are given in the upper part of fig. 1. We 
see that the sum D<^{Q) is in reasonable agreement with pci>{Q) from the Fourier transfor- 
mation of p$(r) inserted in eq. ([3]) (dot-dashed line), which is quite well approximated by 
a radial dependence V'$(r) = tpo exp[~{r /a)'^] with values of k of about 1.5. The resulting 
density is given in the lower part of fig. 1, which indicates clearly that a self-stabilized 
2-gluon field is generated. The mass m$ in the relation between Q and Q' has been used 
as a fit parameter; for a 2-gluon system with a mean square radius of about 0.5 fm^ this 
yields m<j, ~ 0.68 GeV . This is consistent with the gluon pole mass of 0.64 ± 0.14 GeV 
deduced in ref. p]. We shall see, that the extracted mass can be understood as binding 
energy of the 2-gluon system including relativistic mass corrections. 

The extracted 2-gluon density should give a significant contribution to the gluon 2-point 
functions extracted from lattice QCD simulations in form of gluon field correlators [ini [H] 
and the QCD gluon propagator (see e.g. Il2j). In the upper part of fig. 2 we make 
a comparison of our results with the 2-gluon field correlator C±{r) of Di Giacomo et 
al. [11], in the lower part with the gluon propagator of Bowman et al. [T3] . 

We get a good agreement with the lattice data (note that in fig. 2 the gluon propagator is 
multiplied by Q^), but we have to add a second 2-gluon component of smaller size, given 
by the dot-dashed lines. 

3. Bindung potential of gluons and 0++ glueball states 

A finite 2-gluon density as shown in fig. 1 may be interpreted as a bound state of the two 
gluons (glueball). Therefore, from the 2-gluon density the binding potential of the 2-gluon 




(4) 
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Figure 2: Gluon field correlator log{C±) from ref. [TT] (upper part) and gluon propagator 
from ref. [I3] (lower part) from lattice QCD simulations in comparison with our results. 
The lower solid lines correspond to the density p$(r) and its Fourier transform p<s>{Q), 
respectively, and the dot-dashed lines to an additional vector component. The sum of 
both contributions yields a good description of both lattice data. 
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radius (fm) 

Figure 3: Upper part: deduced 2-gluon binding potentials (dot-dashed and dashed hues) 
in comparison with the confinement potential from lattice gauge calculations [1] (upper 
part). The solid line corresponds to the binding potential for the nucleon discussed below. 
Lower part: potentials deduced for heavy vector {ss, cc, and bb) systems discussed below 
in comparison with lattice QCD results [3]. 
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system can be obtained by solving a three-dimensional reduction of the Bethe-Salpeter 
equation in form of a relativistic Schrodinger equation 

where ^^^(r) is the 2-gluon wave function and /i$ a relativistic mass parameter, which is 
related to m$ by yU<j, = ^ m$ + 5m, where 6m is a relativistic correction. Slightly different 
solutions of the binding potential were obtained, which are given by the dot-dashed and 
dashed lines in the upper part of fig. 3. Because of the relation 2g {qq)"' this potential 
can also be considered as confinement potential between the emitted quarks. This is in 
surprising agreement with the 1/r + linear form expected from potential models [H E] 
and consistent with the confinement potential from lattice QCD [1]. It is important to 
note, that our potential reproduces the 1/r + linear form without any assumption on its 
distance behavior; this is entirely a consequence of the deduced radial form of the 2-gluon 
wave function. 

Bound state energies Ei=0.68±0.10 GeV, 1.70±0.15 GeV, and 2.58±0.20 GeV have been 
extracted. Further, in the q-q potential ([3]) we find one bound state with an energy in the 
order of -10 MeV. This very low binding energy indicates clearly, that the glueball states 
must have a large width, since we expect F ~ I/Eq. This is consistent with the general 
expectation for the width of glueball states (F > 500 MeV). From these results we may 
conclude, that the glueball ground state with i?o=0.68±0.10 GeV and a large width may 
be identified with the scalar o"(600). 

Glueball masses have been deduced also from lattice simulations [31 [Tl], in which a glueball 
mass below 1 GeV has not been found. However, in these simulations 0"^^ glueball masses 
have been extracted at about 1.7 and 2.6 GeV, which correspond very nicely to the first 
and second radial excitation in table 1. Our evidence for a low lying glueball is supported 
by QCD sum rule estimates [15], which also require the existence of a low lying gluonium 
state below 1 GeV. 

2. Heavy flavour neutral systems and as{Q) 

Self- consistent 2-gluon densities have been deduced also for smaller systems corresponding 
to the size of ss, cc, and bb mesons. Resulting momentum distributions and the corre- 
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spending Fourier transformed densities are given in fig. 4, which are in good agreement. 




Figure 4: Momentum distributions (histograms) and Fourier transforms of the 2-gluon 
density (dot-dashed hues) for 4 different densities investigated, corresponding to hght qq 
system (a), ss (b), cc (c), and bb (d). The dot-dashed histograms are simulations assuming 
c and b quark masses of 1.4 and 4.5 GeV, respectively. 

It is interesting to investigate the effect of quark masses in our simulations. Using quark 
masses of 1.24, and 4.5 GeV for c and b quarks, respectively, yields the lower dot-dashed 
histograms, indicating that self-consistent solutions are not possible. Thus, for all systems 
the intrinsic quark masses have to be zero (or very small). The 2-gluon binding potential 
is given in the lower part of fig. 3 in good agreement with the confinement potential from 
lattice QCD [% 
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Since all intrinsic quark masses have to be small in our approach, the masses of the 
different systems have to be explained in a different way. Whereas the binding potential 
(of 2-gluons) gives rise to binding energies in the order of 1 GeV, the binding potential of 
quarks ([T]) depends strongly on the size of the 2-gluon densities. Therefore, the binding of 
quarks can be much larger. This is shown for the different systems in fig. 5. For the heavy 
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Figure 5: Folding potential ([T]) for the four different cases (a) - (d) corresponding to light 
glueballs, ss, cc, and bb with the binding energies indicated. 

systems (which are of small size) the binding energies are in the order of 2.4 and 9.0 GeV, 
respectively, which shows that indeed the masses of all systems can be explained by the 
binding of quarks and gluons. The resulting energies of ground and radial $, ^ and T 
states are in good agreement with the experimental spectra. Details of these calculations 
will be given elsewhere. 

The Fourier transformed potential ([2]) is directly related to the strong coupling as{Q)- 
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Using the different 2-gluon density distributions Pi{Q) deduced from fig. 4 we obtain 
as = Sj aip^.{Q). This gives a quantitative description of up to large momenta, see 
fig- 6. 



Strong coupling a,(Q) 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

Q (GeV/c) 



Figure 6: Strong couphng from lattice QCD [16] and experiment |T7j (triangles) in 
comparison with our results, given by the solid line. The contributions of the different 
2-gluon densities pq,XQ) s-^e given also. 

5. Nucleon structure 

Baryons may be described in our approach assuming the decay Ag —>■ 5{qq) —>■ {3q qq) + 
(3g gg), which means 4g — > (baryon + antibaryon). Thus, we decribe the nucleon by 3 
valence quarks coupled to a 2-gluon field yielding pAr(r) = 47r / p3q(r')p$(r — r')dr' . The 
resulting binding potential is given in the upper part of fig. 3 by the solid line, which is 
more shallow than the confinement potential, but the attraction between the emerging 
quarks is increased by a factor 9. The binding energies of the nucleon g.s. and radial 
excitations are 0.94 GeV, 1.42±0.07 GeV, and 1.82±0.12 GeV in good agreement with 
experiment. 
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5.1. Compressibility: from nucleon to nuclear matter 



Finally, we discuss the nucleon compressibility, which may be linked to that of nuclear 
matter. From the excitation of the first radial state, the "breathing mode" of the nucleon, 
the compressibility K]<q has been extracted by operator sum rules |T8l |19] yielding values 
of about 1.3 GeV. The breathing mode has also been investigated in high energy p-p and 
TT-p scattering |20j. From a comparison of transition densities deduced from inelastic p-p 
and e-p scattering strong multi-gluon contributions were extracted, which were about a 
factor 4 stronger than those of the valence quarks. This is in good agreement with the 
present results. From the multi-gluon potentials deduced in ref. [20] we may derive the 
nucleon compressibility directly. Calculating a potential density for the nucleon given by 
y PN{.f) = \ Vnn J PN{rj\f)tj^]\f{r — r^)dr]^ and adding a kinetic energy term T we obtain 
the energy density Epj^ir) = —{V + T) ■ Pn{j). The compressibility is then given by 



KN=r \r=r.. (6) 



From the analysis of high energy p-p scattering [20] the multi-gluon potential is well 
determined. So, we can determine the compressibility. This is given on the left side of 
fig. 7. Indeed, we obtain a compressibility of about 1.3 GeV consistent with the value 
obtained from sum rules [T9l. 



For the case of nuclear matter we assume that the dominant contribution is due to com- 
pressibility of the nucleons (the compressibility due to the binding of nucleons should not 
be much larger than their binding energy). Then the compressibility is related to the 
central (scalar) nucleon potential, which has a mean square radius >1.5 fm^. Using the 
corresponding density the derived compressibility is in the order of 160-170 MeV, this 
is shown on the right side of fig. 7. This value is rather close to the compressibility of 
nuclear matter Kao of about 220-250 MeV deduced from the study of the giant monopole 
resonance in heavy nuclei. 

6. Conclusion 



The present solution of the confinement problem based on our phenomenological descrip- 
tion of two-gluon fields differs entirely from earlier suggestions, that confinement could 
arise from complicated non-perturbative field configurations (magnetic monopoles, fiux 
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Figure 7: Left: Nucleoli potential density (solid line) and derived compressibility function 
([6]) (dot-dashed line). Right: The same for the scalar nucleon potential, important for 
nuclear matter. 

tubes, vortices or strings) in the Abelian projection of QCD, which had severe problems, 
e.g. with Casimir scaling. Further, none of these models could explain the generation of 
mass and the complex Yang-Mills gluon structure found in lattice QCD calculations. In 
our description all these problems are tied together and are well described. The masses 
of hadrons are explained by binding effects, and all intrinsic quark masses have to be 
consistent with zero. Thus, a scalar Higgs field in which the quark masses are generated 
is not needed. Also the axion problem does not exist when quark masses are zero. 

It is of large interest that in our approach, in which many of the properties of hadrons are 
well described, including the problem of the light pion mass and the non-existence of chiral 
symmetry, quarks arise only from the decay of multi-gluon systems. This has important 
consequences for our understanding of the origin of our universe and of baryogenesis. 

Finally, by the discussion of the compressibility a first example is given, which shows that 
the properties of hadrons are strongly tied to those of nuclear systems. 
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